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1.  INTRODUCTION 

In  antiproton  annihilation  propulsion  of  spacecraft,  a  possible  choice 
for'the  annihilation  process  is  low-energy  antiprotons  (kinetic  energy  <  1  eV) 
annihilating  with  the  protons  of  hydrogen  atoms.^  At  such  energies,  an 
atomic  rearrangement  occurs  in  which  the  electron  (e-)  of  the  hydrogen  atom 
(H)  is  emitted  while  the  proton  (p,  positive)  of  the  hydrogen  atom  and  the 
antiproton  (p  ,  negative)  form  an  excited  bound  state  of  protonium  (Pn)^: 

p  +  H  -*■  Pn  «■  e“  ( 1 ) 

The  p  and  p  in  Pn  inevitably  annihilate  after  undergoing  radiative  decay  to 
states  of  lower  energy.  The  energy  lost  by  the  p  by, entering  a  negative- 
energy,  bound  state  is  transferred  to  the  e7  That  energy  is  slightly  in 
excess  of  the  13.6  eV  binding  energy  of  the  e“in  H.  The  cross  section  for 
annihilation  through  reaction  1  at  energies  below  about  20  eV  is  much  greater 
than  that  for  annihilation  of  an  antiproton  on  a  bare  proton. Hence  the 
interest  in  this  and  similar  reactions  for  U3e  in  antiproton  annihilation 
propulsion.  Similarly  high  annihilation  cross  sections  result  when  the  H  is 
replaced  by  another  atom  or  molecule  and/or  the  p  is  replaced  by  an 
antihydrogen  atom  (an  antiproton  with  a  positron  bound  to  it,  which  is  emitted 
along  with  the  electron). (2>3) 

1.1  Annihilation  Cross  Section 

Morgan  and  Hughes  determined  the  cross  section  for  reaction  1  (see" 

Fig.  1)  by  employing  a  semiclassical ,  impact-parameter  approximation  in  which 
the  p  -  H  relative  motion  is  assumed  to  be  classical,  while  the  motion  of  the 
e“  in  H  is  treated  quantum  mechanically.^2’^  The  e"  motion  is  determined  by 
solving  the  Schroedinger  equation  for  the  ground  state  of  the  e“  in  the  field 
of  a  negative  charge  (the  p  )  and  a  positive  charge  (the  p)  that  are  separated 
by  a  fixed  distance,  R  (see  Fig.  1).  The  approximation  of  holding  the  p 
and  p  (or  other  nuclei)  fixed,  while  determining  the  electron  motion,  is  known 
as  the  fixed-nucleus  or  Born-Oppenheimer  approximation.^  Solving  for  the 
electron  motion  gives  the  electron  energy  as  a  function  of  R.  This  energy, 
along  with  the  p  -  p  attraction  determines  the  potential  energy  (or  force) 
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Fig.  1.  Morgan-Hughea  description  of  p  -  H  scattering  (not  to  scale). 

a:  R  is  held  fixed  while  the  e~  motion  is  determined  quantum 

mechanically.  This  provides  the  potential  energy  between  the 
p  and  the  H  atom  to  determine  the  classical  orbits  of  the  p  . 

b:  When  the  impact  parameter,  b,  of  the  p  i3  >  bQ  (bc  depends  on  the 
p  initial  kinetic  energy)  then  Rmin  >>  Rc. 

c:  When  b  <  bc  then  Rmin  <<  Rc  and  rearrangement  occurs  with  the  e~ 
emitted  and  the  p  becoming  bound  to  the  p.  The  annihilation 
of  p  and  p  then  follows. 

i 


Parts  b  and  c  are  on  the  following  page. 
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between  the  p  and  the  H,  which  is  used  in  calculating  the  relevant  features  of 

the  p  -  H  classical  orbits  (R  not  fixed). (2.3) 

The  classical  orbits  have  an  interesting  property  for  antiproton  kinetic 

energies  less  than  about  1  eV.  When  the  inpact  parameter,  b,  is  greater  than 

a  particular  value,  bc,  which  depends  on  the  collision  energy,  the  inner 

turning  point  (point  of  closest  approach,  designated  Rm^n)  is  not  much  smaller 

than  bc;  the  p  and  H  stay  relatively  far  apart  and  no  rearrangement  (reaction 

1)  occurs  (see  Fig.  1).  However,  as  b  becomes  less  than  bc,  Rrain  drops 

discontinuously  from  large  values  to  be  much  smaller  than  the  critical  radius, 

Rc.  Rc  is  the  value  of  R  less  than  which  the  e  is  no  longer  bound  to  the  p 

in  H  because  of  the  strong  repulsive  force  of  the  p  .  Wightman^  determined 

the  value  of  RQ  to  be  0.639aQ  where  aQ  =  5.29  x  10-^  cm  is  the  first  Bohr 

radius  of  hydrogen.  (Wightman's  calculation  was  for  negative  muon  -  hydrogen 

which  is  equivalent  to  p  -  H  in  the  Born-Oppenheimen  approximation). 

Following  Wightman,  Morgan  and  Hughes  assumed  that  once  R  became  less  than 

Rc,  rearrangement  occurred  with  near-unity  probability;  the  electron  left  the 

region  and  the  p  became  bound  to  the  p  since  it  had  lost  energy  to  the 

escaping  electron.  Thus  the  rearrangement  and  annihilation  cross  section  was 
2 

o  =  irb  .  This  cross  section  is  shown  in  Fig.  2. 
c 


Morgan  and  Hughes  demonstrated  the  validity  of  using  classical  orbit3  for 

the  p  -  H  motion  by  showing  that  the  results  were  not  altered  when  that  motion 

(  ?  ) 

was  treated  quantum  mechanically  with  an  optical-model  potential  energy.  ' 
However,  their  arguments  that  with  high  probability  the  e"  left  the  scattering 
region  permanently  for  b  <  bc  and  did  not,  in  fact,  reattach  itself  to  the 
proton  with  a  significant  probability  were  only  qualitative.  If  there  is  a 
possibility  for  reattachment ,  then  the  rearrangement-annihilation  cross 
section  is  lower  than  their  value.  The  problem  in  considering  reattachment  is 
that  the  Born-Oppenheimer  approximation  breaks  down  for  R  £  Rc.  Thus,  in  lieu 
of  a  full  quantum-mechanical  solution  of  the  problem,  with  essentially  no 
approximations ,  it  is  difficult  to  quantitatively  demonstrate  that  the  e“  is 
indeed  emitted  with  high  probability  when  R  becomes  less  than  Rc.  The  purpose 
of  the  work  reported  here  has  been  to  examine,  quantitatively,  what  goes  on 
during  the  time  when  the  Born-Oppenheimer  approximation  breaks  down. 


/ 
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2.  THE  p  -  H  SYSTEM  AND  BORN-OPPENHEIMER  BREAKDOWN 

The  Born-Oppenheimer  approximation  is  valid  when  the  mean  velocity  of  the 
electron  relative  to  the  proton  is  much  larger  than  the  velocity  of  the 
antiproton  relative  to  the  proton.  Such  is  true  at  the  kinetic  energies  of 
interest  when  R  is  a  few  times  aQ  or  greater,  principally  because  the 
antiproton  mass  (same  as  the  proton  mass)  is  about  1836  times  the  electron 
mass.  However,  as  R  decreases,  the  mean  velocity  of  the  electron  decreases  as 
it  is  pushed  farther  from  the  proton  and  its  energy  approaches  zero  from  a 
negative  value  while  the  velocity  of  the  antiproton  increases  due  to  its 
attraction  to  the  H  atom.  For  R  =  RQ  in  the  Born-Oppenheimer  approximation, 
the  energy  and  velocity  of  the  electron  have  dropped  to  zero,  so  the  cross¬ 
over  point  of  the  velocities  is  at  a  small  value  of  R  but  greater  than  RQ. 

I  ' 

Before  considering  what  happens  when  the  Born-Oppenheimer  approximation 

breaks  down,  it  is  necessary  to  determine  the  value  of  R,  in  that 

approximation,  at  which  the  cross-over  point  of  the  velocities  occurs.  That 

value  of  R,  designated  Rx,  was  found  by  using  two  different  methods  to 

determine  u,  the  mean  velocity  of  the  electron.  Both  make  use  of  the  Born- 

Oppenheimer  energy  of  the  electron,  e,  which  is  given  in  Ref.  2,  (add  aQ/R  -  V2 

to  V  in  Table  VI  of  that  reference  to  obtain  e  in  atomic  units),  where  for 

( S ) 

values  of  R  less  than  5  aQ,  e  is  taken  from  work  by  Wightman  and 
Bates. ^  e  is  given  in  Table  1  along  with  V,  the  p  -  H  potential  energy. 
e  and  V  are  related  by 

V  -  e  -  eVR  -  eQ  ,  (2) 

where  e  =  4.80  x  10  10  cm^/2  g1/2  s  1  is  the  unit  electric  charge,  and 

e  »  -V2eJ/a  is  the  ground  state  energy  of  the  electron  in  H.  In  the  first 

(7) 

method,  the  virial  theorem  is  assumed  to  apply  to  e,  so  the  kinetic  energy 
of  the  electron  is  -  V2  e  and  u  is  therefore  (-e/m)  2  ,  where  m  =  9.11  x  10  g 
is  the  electron  mass. 


Table  1.  Various  enemies  pertinent  to  the  p  -  H  system.  t.  an.l  V  are  taken 
from  References  2,  5,  and  6.  aQ  -  one  atomic  unit  of  distance  = 

*>  5.29  x  10  9  cm,  e2/a  =  one  atomic  unit  of  energy  =  27.2  eV. 


p  -  p 

separation, 

RCaJ 


P  ~  H 

potential 

energy, 

V[e2/a?] 


e  energy, 
e[e2/a_]  , 


adj  us ted 
e“  energy, 
e '  [  e  2  /  a  .  ] 


-9.500 

0 

.■'1 

0  ? 

-4.500 

0 

0  i 

-2.000 

0 

0  i 

S' 

-1.565 

0 

-0.00038* 

0  ! 

0.0043*  | 

1 

% 

-0.929 

§ 

-0.510 

-0.010 

-0.183  < 

8 

-0.037 

1 

-0.037 

i 

-0.412  J 

-0.0091 

-0.2591 

-0.4909  1 

-0.00136 

-0.3580 

-0.500013  J 

§ 

-0.000310 

-0.400310 

-0.500183  i 

i 

-0.0000141 

-0.4500141 

-0.5000140 

t 

1 

*  interpolated  using  e  or  e'  »  const, 


x  (R/a0  -  0. 639) v,  v  =  2 


A  problem  with  determining  u  in  the  first  method  involves  the  fact  that 
e  has  a  substantial  contribution  from  the  potential  energy  between  the  p  and 
the  quantum  mechanically  distributed  e~.  For  values  of  R  equal  to  a  few  aQ 
(and  greater)  this  is  mainly  a  static  contribution;  it  has  only  a  small  ^effect 
on  the  dynamics  of  the  electron  motion.  Hence  the  second  method  for 
determining  u,  in  which  part  of  the  p  -  e“  potential  energy  is  subtracted 
from  c  (making  it  more  negative)  before  applying  the  virial  theorem.  The 
fraction  of  this  potential  energy  to  subtract  is  uncertain,  but  the  fraction 
is  about  unity  when  R  is  large  and  less  than  unity  for  R  =■  1  or  2  aQ  where 


the  p  affects  the  dynamics  of  the  electron  motion.  For  the  fraction, 

-R/a 

1  -  e  is  chosen  where  R  is  in  units  of  aQ  and  a  is  the  mean  radius  of  the 
e"  distribution  (also  in  units  of  aQ).  Thus, 


where  the  first  quantity  in  paranthesis  is  the  fraction,  ai  .i  he  second 
quantity  in  parenthesis  (with  an  inner  parenthetic  expression)  is  the  mean 
p  -  e“  potential  energy.  In  Eq.  3,  R  is  in  units  of  aQ  and  e  and  e1  are  in 
units  of  eJ/a0.  The  p  -  e“  mean  potential  energy  was  determined  by  assuming 

that  the  e"  distribution  is  approximately  hydrogenic,  i.e.  it  has  the  form 

const,  x  e-3R/a  (wave  function  =  const,  x  e~3R/(2a)^  Under  the  same 

V, 

assumption,  a  =  (-9/(8  e'))  -.  When  this  equation  for  a  is  substituted  into 

Eq.  3,  that  equation  can  be  solved  iteratively  for  e’  as  a  function  of  R.  The 

resultant  values  of  e’  are  given  in  Table  1. 

The  values  of  c'  for  R  >  2  aQ,  where  the  p  is  outside  of  the  e“ 
distribution,  appear  to  represent  the  dynamic  motion  of  the  e~  more  reliably 
than  the  values  of  e  since  they  are  closer  to  the  value  of  -V2Which  both 
e  and  e1  have  for  R  =  ®.  For  smaller  R  than  2  aQ  it  is  surprising  that  s’  is 
so  much  greater  in  magnitude  than  e.  Thus  it  is  likely  that  a  better  value 
of  e’  probably  lies  between  e  and  the  present  value.  As  will  be  seen  below, 
this  uncertainty  in  e'  has  only  a  small  effect  on  Rx  because  of  the  strong 
dependence  of  u  on  R, 

Applying  the  virial  theorem  to  both  e  and  e'  yields  the  values  for  the 
mean  electron  velocity,  u  -  (-2e/m)T/2  ,  or  the  same  with  e'  in  place  of  e, 
where  m  is  the  mass  of  the  electron,  that  are  given  in  Table  2.  Also  given  in 
Table  2  i3  the  velocity  of  the  antiproton  relative  to  the  proton.  This  latter 
quantity  is  determined  for  an  antiproton  collision  energy  of  0.001  ez/aQ  » 
0.027  eV.  It  applies  for  all  impact  parameters  that  are  less  than  bc.  The 
results  are  essentially  unaltered  near  R  -  R  if  lower  or  higher  collision 
energies  are  used  a3  long  as  the  energy  is  <  1  eV.  (At  higher  energies  the 
method  used  by  Morgan  and  Hughes^2’^  to  determine  the  cross  section  may  not 
apply  because  the  discontinuity  of  Rmin  becomes  considerably  less  pronounced 
or  disappears). 

By  plotting  the  velocities  in  Table  2,  it  is  found  that  Rx  -  0.79  aQ  when 
e  is  used  and  Rx  -  0.68  aQ  when  e'  is  used  -  nearly  equal  values  in  spite  of 
the  large  differences  between  e  and  e'.  Since  the  Born-Oppenheimer 
approximation  becomes  invalid  as  R  decreases  through  Rx,  it  will  be  assumed 
that  it  is  valid  until  R  »  Rx.  In  the  following  section,  the 
evolution  of  the  p  -  H  system  is  modeled  for  R  <  Rx. 


Table  2.  Values  of  the  electron  mean  velocity,  u  and  u ' ,  compared  to  the 
antiproton  velocity,  v.  u  and  u'  are  obtained  by  applying  the 
•>  virial  theorem  to  the  electron  energy  e  and  the  adjusted  electron 
energy  e',  and  v  is  for  an  initial  collision  energy  of  0.027  eV. 
x  Units  for  the  velocities  are  e(aQm)  ^(atomic  velocity  units). 


RCaJ 

u 

u' 

V 

0. 1 

0 

0 

0.2012 

0.2 

0 

0 

0..  1  385 

0.4 

0 

0 

0.0924 

0.639 

0 

0 

0.0674 

0.7 

0.028 

0.092 

0.0630 

1  .0 

0.1  42 

0.605 

0.0467 

2.0 

0.417 

0.907 

0.0194 

4.0 

0.720 

0.991 

0.0066 

7.0 

0.847 

1.000 

0.0032 

10.0 

0.895 

1.000 

0.0024 

20.0 

0.949 

1.000 

0.0021 

3.  THE  p  -  H  SYSTEM  FOR  R  £  Rx 

For  R  i  Rx  the  value  of  e  in  Eq.  2  for  the  p  -  H  potential  energy  V  is 
small  enough  to  be  neglected,  so  for  S  S  Rx 


e2  q  2 

V  a  -  _  +  S - 

R  2a 


where  eQ  has  been  replaced  by  its  value,  -V2e2/a0.  Because  t  (or  e'  )  is  so 
small,  the  mean  radius  of  the  electron  distribution  is  much  greater  than  Rx. 
Thus,  V  is  more  properly  termed  the  p  -  p  effective  potential  for  R  £  R^.  It 
is  equal  (c.f.  Eq.  4)  to  the  p  -  p  potential  energy,  -e2/R  plus  the  energy, 
V2e2/aQ,  that  has  been  transferred  from  the  antiproton  to  the  electron.  This 
energy,  equal  to  13.6  eV,  is  the  negative  of  eq.  The  initial  kinetic  energy 
of  the  antiproton  (£  1  eV)  is  less  than  the  energy  transferred,  hence  the  fact 
that  the  p  will  become  bound  to  the  p  (in  a  negative  energy  state)  unless  the 
e"  rejoins  the  p  and  transfers  the  absorbed  energy  back  to  the  p  .  Due  to  the 
simple  form  of  Eq.  4  and  because  it  involves  tire  inverse  of  R,  analytic 
formulae  for  the  orbital  motion  of  the  p  relative  to  the  p  can  be  determined 
for  R  <  Rx.  In  addition  approximate  means  exist  for  following  the  motion  of 
the  e"  during  this  time  period. 
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The  analytic  formulae  for  the  motion  of  the  p  are  the  standard  formulae 
for  an  attractive  1 /R  potential  (1/R2  force)  but  with  an  energy  ofV2e2/aQ 
subtracted  from  the  initial  p  kinetic  energy.  These  formulae  give 


t 

x 


(2M) 

1-2K 


2  J/. 

-  ( ( 1 “2K  )  ' 


'2(1 


.  1  -  R  (1-2K  ) 

sin  ( - ij) ) 

( 1 -2b2K  (1-2K  ))  2 
o  o 


-(-R2  (1-2K  )  *  2R 
x  o  X 


2b2K  ) 
o 


)h 


(5) 


•VJ 


\ 


and 


R  -  2b2K 

9  =  ir  +  2  sin  (— - - - - jj) 

x  R  ( 1 -2b2 ( 1 -2K  ))  2 

x  o 


(6) 


for  the  time  that  R  is  %  Rx  and  for  the  angle  (centered  on  the  p)  through 
which  the  p  moves,  while  R  2  Rx,  in  its  orbit  around  the  p.  In  Eqs.  5  and  6, 

Kq  is  the  initial  collision  kinetic  energy  of  the  p  (  <  1  eV)  in  the 

laboratory  frame,  i.e.,  the  p  has  no  initial  kinetic  energy.  (In  the  center 
of  mass  frame  the  initial  kinetic  energy  is  KQ/2).  In  Eqs.  5  and  6,  M  -  1836 
i3  the  p  (»p)  mas3  in  units  of  m  (the  electron  mass),  Rx  is  in  units  of  aQ,  KQ 
is  in  ’units  of  e2/aQ,  tx  is  in  units  of  aQ^2  e~'  (atomic  time  units), 
and  9X  is  in  radians.  Figure  3  depicts  the  geometry  of  the  p  -  p  motfon  for 
R  £  R  .  For  K  *  0  to  1  eV  and  b  *  0  to  b.  values  of  t  range  from  23  to  27 
atomic  ’units  for  Rx  »  0.79  aQ  and  from  18  to  22  atomic  units  for  Rx  -  0.68  aQ, 
with  the  maximum  values  of  tx  being  for  the  maximum  values  of  KQ  and  b,  and 
the  minima  of  tx  being  for  KQ  »  0  and  b  »  0.  For  the  same  ranges  of  KQ  and 
b,  0X  varies  from  303°  to  360°  for  Rx  -  0.79  aQ  and  from  296°  to  360°  for  Rx  - 

0.68  aQ,  with  the  maximum  values  of  9^  (360°,  which  means  back  to  the  same 

point  after  a  full  circuit)  coming  from  b  -  0  (any  KQ  in  the  range)  and  the 
minima  of  9X  from  the  maximum  values  of  KQ  and  b. 
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The  e-  motion  is  approximated  as  that  of  an  evolving  Gaussian  wave 
packet,  during  the  time  that  R  S  Rx,  while  the  p  -  p  configuration  is  changing 
as  described  above.  It  is  assumed  that  the  evolution  is  unaffected  by 
the  p  and  p  because  of  the  large  spread  of  the  e~  distribution  compared  to  R 
for  R  £  Rx.  Thus,  the  wave  function  for  the  e~  is 


V  ~  -  m“3/4(3  +  it/S)“3/2  e"  2( 82  +  it ) 


where  t  is  the  time  following  the  first  instance  that  R  -  Rv,  B  is  a  constant, 

1  * 

all  quantities  are  in  atomic  units,  and  i  -  (-1r2.  Equation  7  is  the 
solution  to  the  time-dependent  Schroedinger  equation 


-V  s  +  * 


The  electron  probability  distribution  (referred  to  above  as  just  the  electron 
distribution)  is  the  square  of  absolute  value  of  the  wave  f 'unction, 


I'fg-I2  -  iT3/2(e2  ♦  t2/B2)“3/2  e“  B2+t2/B; 


The  kinetic  energy  of  the  electron  in  this  description  is  3/ C ^S2 ) .  This 
energy  may  be  equated  to  the  kinetic  energy  that  the  electron  ha3  at  t  -  0 
(R  =*  Rx)  (when  the  Born-Oppenheimer  approximation  breaks  down).  Using  the 
virial  theorem,  this  latter  energy  is  -e/2  or  -e’/2  at  R  -  Rx.  (For  R  -  Rx, 
e  -  -0.00166  and  e’  -  -0.00211).  Thus 


S  -  (-3/(2e))/2  or  (>3/(2e * )  )'2  , 


and  substituting  the  values  of  e  and  e'  gives  B  -  30.1  aQ  for  e  and  S  -  26.7 
aQ  for  e'.  Hence  the  earlier  statements  that  the  e~  distribution  is  spread 
over  distances  much  greater  than  RY. 


3.1  Reattachment  Probability 

j 

The  likelihood  that  the  e“  will  reattach  itself  to  the  p  once  R  again 
becomes  greater  than  Rx  is  related  to  the  similarity  between  two  forms  of 
i|»e-  .  The  first  form  is  ijie-  at  t  »  t  ,  when  the  Born-Oppenheimer 
approximation  can  be  re-established  and  the  e~  motion  can  again  become  coupled 
to  the  p  -  p  motion,  and  the  second  is  the  \^e-  that  the  electron  would  have 
if  it  were  to  become  reattached  to  the  p.  From  the  viewpoint  in  the  center  of 
mass  frame  of  the  p  and  p,  the  former  i|>  -  will  be  taken  as  given  by  Eq.  7  for 
t  =  tx,  located  and  moving  according  to  the  overall  velocity  it  had  at  t  *  0, 
while  the  latter  will  be  a  i|»e-  as  given  by  Eq.  7  at  t  =  0  (last  time  of 
attachment)  but  with  the  position  and  velocity  of  the  p  at  t  -  tx.  The 
changes  in  position  and  velocity  may  be  seen  in  Fig. » i4- where  the  p  -  p  motion 
for  R  £  Rx  is  shown  in  the  center  of  mass  frame.  In  the  above  "position" 
means  the  point  in  space  to  which  r  >  0  in  Eqs.  7  and  9  refers,  while  "overall 

velocity"  refers  to  the  velocity  of  that  point  in  the  center  of  mass  frame. 

The  meaning  of  "similarity"  will  be  semiquantitati ve  for  size  and  position 
changes  of  but  will  be  given  a  precise  definition  when  the  effects  of 

velocity-change  are  considered. 

When  the  values  of  B  along  with  t  -  0  and  tx  are  placed  in  Eq.  7,  it  is 

seen  that  there  i3  essentially  no  increase  in  size  or  other  change  in  form  in 

\jie-  due  to  time  evolution  from  t  =*  0  to  t  -  tv.  Since  the  size  of  <|ie-  (the 
mean  radius  is  2ir  2  (B1  +  t2/B2)  2)  is  much  greater  than  the  changes  in 
position  (which  are  not  larger  than  Rx),  changes  in  position  will  not  result 
in  significantly  dissimilar  wave  functions.  Thus,  the  only  possible 
significant  lack  of  similarity  between  the  wave  functions  to  be  compared  is 
that  which  can  come  from  the  difference  in  velocity  of  their  centers.  The 
difference  in  velocity  of  the  centers  of  the  e“  wave  functions  is  the 
difference  between  the  velocities  of  the  proton  in  the  center  of  mas3  system 
at  t  *  0  and  t  »  tx.  Because  9X  is  about  360°  (change  in  direction  about 
180°)  the  change  in  velocity  of  the  proton  in  phe  center  of  mass  system  is 
about  twice  its  velocity  at  R  -  Rx  (t  -  0).  This  velocity  is  one-half  the  p  - 
p  (or  p  -  H)  relative  velocity,  so  the  change  in  velocity  of  the  center  of  the 
electron  wave  function  is  equal  to  v  at  R  -  Rx.  Thus,  the  velocity  difference 
is  equal,  in  atomic  units,  to  0.058  (using  c)  or  0.065  (  jsing  e ' ) . 


Fig.  4.  Motion  of  the  p  and  p  in  the  center  of  mass  frame  for  R  S  R  (not  to 
scale).  At  t  -  0,  when  the  e“  motion  becomes  uncoupled  from  the  p 
motion,  the  e~  wave  function,  'J'e-*  is  centered  at  point  A  (for  the  p) 
and  has  an  overall  velocity  equal  to  the  p  velocity  at  t  *  0.  When 
the  p  has  reached  point  B  (R  again  «  Rx)  the  center  of  i|ie-  has  moved 
to  point  C.  The  reattachment  e-  wave  function  has  a  form  equivalent 
to  that  at  point  A,  but  it  is  centered  at  point  B  and  has  an  overall 
velocity  equal  to  the  p  velocity  at  that  point.  In  the  center  of 
mass  reference  frame,  the  p  and  p  velocities  are  equal  in  magnitude 
and  opposite  in  direction.  The  magnitude  is  half  that  of  the 
relative  velocity  of  the  p  and  p. 
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Because  of  the  small  effects  of  change  in  form  and  position  of  e“,  it  may 
be  taken  to  be 


-3/4  -3/2  • 

4»  -  =  it  3  e 


r' 

2SZ 


(11) 


when  t  =  t  and  at  the  same  time  the  reattached  \J>e-  may  be  described  by 


♦  -  -  S3'*  B~3/2  e"  28 
e 


ikz 


(12) 


where  the  multiplicative  factor,  e  ,  is  due  to  the  change  in  velocity.  The 
quantity  k  is  the  wave  number  for  this  change  and  is  equal  to  v  in  atomic 
units  (since  the  mass  of  the  e~  is  'unity  in  atomic  units)  while  z  is  the 
cartesian  coordinate  in  the  direction  of  the  velocity  change  (the  results 
below  are  not  affected  by  the  choice  of  direction  for  the  z-axis). 

The  similarity  between  the  two  forms  of  4>e-  given  by  Eqs.  11  and  12  will 
now  be  defined  specifically  as  the  square  of  the  absolute  value  of  the  overlap 
integral  of  these  two  forms  for  This  quantity,  in  the  "sudden" 

approximation,  is  the  probability  that  the  e"  will  reattach  to  the  p  as  R 

passes  through  Rx  on  the  way  out.  The  sudden  approximation  is  accurate  when  a 

system  encounters  a  change  in  its  potential  energy  that  occurs  during  a  time 

interval  that  is  much  shorter  than  the  time  required  for  that  system  to ‘change 
significantly.  Such  conditions  are  met  fairly  well  here.  The  result  of  the 
integration  gives 

B2v2 

P  -  e~  2  (13) 

for  the  probability  of  reattachment .  Using  the  values  of  S  and  v  already 
given,  P  -  0.22  when  either  e  or  c'  is  used  to  obtain  8  and  v. 
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At  face  value,  a  value  for  P  of  about  20?  means  that  the  Morgan-Hughes 
results  for  the  p  -  H  rearrangement-annihilation  cross  section'2, should  all 
be  redyced  by  about  20?.  However,  the  accuracies  of  two  approximations  that 
went  into  this  value  are  unclear.  The  first  approximation  is  the  assumption 
that  the  Born  Oppenheimer  approximation  is  valid  and  accurate  for  all  values 
of  R  greater  than  R  ,  the  point  where  the  electron  mean  velocity  is  equal  to 
the  antiproton  velocity,  and  that  the  electron  motion  is  wholly  uncoupled  from 
the  p  and  p  for  R  <  Rx.  The  second  approximation  is  the  assumption  that  the 
electron  wave  function  is  a  gaussian  wave  packet  for  R  <  R  .  In  reality,  the 
transition  from  Born-Oppenheimer  motion  of  the  electron  to  uncoupled  motion 
occurs  gradually  over  a  range  of  values  of  R,  and  the  form  of  the  wave 
function  for  R  <  Rx  is  more  complicated  than  the  gaus3ian  form  assumed. 

< 

4.  CONCLUSIONS  AND  DISCUSSION 

A  detailed  investigation  has  been  made  of  antiproton  -  hydrogen  atom 
scattering  for  small  values  of  the  antiproton,  hydrogen  atom  separation  where 
the  Born-Oppenheimer  approximation  is  invalid.  The  calculations  yield  an 
estimate  of  20?  for  a  reduction  in  the  rearrangement  cross  section  calculated 
by  Morgan  and  Hughes^2'-^  that  lead3  to  antiproton  annihilation  with  the 
proton  of  the  hydrogen  atom  at  antiproton  energies  of  several  eV  or  less. 

Such  a  reduction  is  of  limited  negative  consequence  at  the  present  time 
to  considering  antiproton  annihilation  as  a  means  of  spacecraft  propulsion. 

It  means,  for  instance,  that  the  density  of  hydrogen  atoms  (or  other  forms  of 
matter  for  which  similar  reductions  would  probably  occur)  in  the  annihilation 
region  would  have  to  be  increased  by  only  20$  over  previous  estimates^  to 
achieve  the  same  annihilation  efficiency.  However,  because  of  the  approximate 
nature  of  the  present  calculations  and  because  precise  values  of  the 
annihilation  cross  section  will  be  required  in  the  future,  it  is  important  to 
recalculate  the  antiproton  -  hydrogen  atom  rearrangement  cross  section  with  an 
accurate,  f ully-quantum  mechanical  method. ^ 
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